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Time correlation functions of a hybrid exotic meson operator, with J^'^ = 1 ^ , generated in 
quenched lattice QCD are subjected to a (Bayesian) maximum entropy analysis. Five distinct 
spectral levels are uncovered. Their extrapolation into the physical pion meiss region suggests a 
possible relationship to experimentally known states 7ri(1400) and 7ri(1600), and also to a state in 
the 2 GeV region carrying the same quantum numbers. 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) in principle permits 
the existence of hadrons containing both quarks and va- 
lence gluons. These are commonly known as hybrids. 
For example, in the meson sector, the presence of explicit 
glue allows for states with quantum numbers J^'" = 1 ^. 
The standard quark model cannot describe a meson with 
the above quantum numbers. Hadrons with this feature 
are called exotic. Other exotic states are possible 0, 
however, in this article we will solely consider the 1 ^ 
hybrid exotic meson. 

According to Q experimentally established 1 ^ 
mesons are 7r(1400) and 7r(1600). Several experimental 
programs currently focus on studies of exotic hadrons in 
general. 

On the theoretical side we expect the lattice formula- 
tion of QCD, being a first-principles approach, to eventu- 
ally yield the properties of hybrid exotics such as masses 
and decay widths. Studies in this area have been but a 
few [3., 4, 5, 6(, |3[J1, 0, ^3 ■ Results hint at masses in the 
1.4 GeV region M, and also at a state close to 1.9 GeV 

In the present article we examine the mass spectrum 
of 1 ^ mesons based on lattice QCD. We utilize the out- 
come of a numerical simulation, gauge field configura- 
tions and quark propagators, which had been previously 
generated to study hybrid exotic decay widths 0. 

The novel feature here is the use of Bayesian analysis, 
more precisely the maximum entropy method (MEM) 
[Tl| . This analysis tool enables us to utilize the entire 
range of lattice time slices, thus eliminating the neces- 
sity for a subjective choice of a time interval defining 
a plateau in the effective mass function, an approach 
almost universally used as an analysis method to ex- 



tract hadron masses. When all time slices are utilized, 
on a given lattice volume, the method delivers unam- 
biguous results, free from subjective choices of effective 
mass plateaus. In addition, the MEM is quite naturally 
adapted to multiple states, if present. We will particu- 
larly rely on this feature to uncover the 1 hybrid exotic 
meson mass spectrum, up to somewhat above 2 GeV. 



II. LATTICE SIMULATION 

In the context of the decay width calculation "t] it was 
argued that systematic errors would likely dominate the 
final result. For this reason the lattice action and lattice 
sizes were selected to keep the numerical effort moder- 
ate. Although such a constraint would not seem to be 
necessary for our present aim of 'just' extracting spectral 
masses, we find that the lattices and quark propagators 
generated in give surprisingly relevant results. 

Thus, we here present the 1 ^ exotic hybrid mass spec- 
trum based on the Wilson gauge field action used with 
Wilson fermions, in the quenched approximation, on an 
anisotropic lattice of size 12^^ x 24. The bare aspect ratio 
is tts/ at = 2, with being the spatial and at the tempo- 
ral lattice constants. Simulations were done at one value 
of a global gauge coupling /3 = 6.15 and four values of 
global hopping parameters they appear in Tab. ^ 

below. We refer the reader to for a precise definition 
of these global parameters. 

In order to set the physical mass scale, standard local 
zero-momentum meson operators of the form 



Ox{t) 



] daixt)TUa{xt) , 



(1) 
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where a is color, were employed for the pseudoscalar me- 
son A = TT, F = 75, the vector meson X = p,r — ji, and 
the oi meson X = ai,F — "fij^. where i — 1,2,3 are 
spatial directions. 

For an operator coupling to 1 ^ hybrid meson states, 
which we collectively refer to with h, we follow Isl and 
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use 

l<i<j<3 X 

(2) 

Here a, b denote color indices and F^j (x) is a product of 
SU{3) link matrices arranged in a clover pattern 

F^.,{x) = U^,{x)U,{x + fL)Ul{x + vWlix) 
+ U„{x)Ul{x + v)Ul{x - fi)Uf,{x - /i) 
+ Ul^{x - jl)Ul{x - jl - i')U^{x - fi- i')Uy{x - v) 
+ Ul{x - u)U^{x - u)U,{x + fi - u)Ul{x) , (3) 

which is used in the spatial planes only thus employing 
magnetic type gluons, in the rest frame. 

The parity transformation V applied to ^ gives 
VOh{t)V~^ = ~Oh{t), as it should. This relation re- 
lies on VUi{x,t)V~^ = U-i{—x,t) for i = 1,2,3, using 
the notation of Ref. ^2 • the quantum level the cor- 
responding propagator respects exact (negative) parity 
because the (quenched) gauge field action S[U] is transla- 
tionally invariant, and satisfies VS[U]V^^ — S[U]. How- 
ever, in a numerical simulation those two properties are 
achieved only in the limit of a large number of gauge field 
configurations, but are otherwise approximate. Conse- 
quently, one should expect that the hybrid meson prop- 
agator is contaminated by states of the wrong (positive) 
parity, albeit with a small amplitude. We will revisit this 
point in the discussion of results. 

With regard to charge conjugation C we encounter a 
similar situation. For the purpose of discussion, chang- 
ing /i to a charge neutral operator /i*^ {du — » dd, uu) , the 
proof of CO;,o(t)C-i = -Oho{t) relies on CF,j{x)C-^ = 
F*j{x), which comes from CUi{x)C^^ = U*{x). Again, 
at the quantum level the corresponding propagator re- 
spects exact (positive) charge conjugation because of 
S[U] = S[U*]. In contrast to the parity case this relation 
is easily enforced in the simulation. The configurations 
[U] and [U*] are equally probable. Thus with each [U] in 
the ensemble of 200 configurations we also include [[/*] 
and compute fermion propagators for both of those. This 
strategy doubles the number of fermion propagators that 
need to be computed to 400. However, charge conjuga- 
tion is now numerically exact, and this also appears to be 
the reason for an observed noise reduction of simulation 
signals. 

All meson operators Ox{t), X ~ tt, p, oj, /i, are ex- 
panded by employing quark field smearing |l3| and gauge 
field fuzzing [l^. 

Smeared quark fields are constructed by spreading the 
original field ipi^) over neighboring sites using a recursive 
procedure. Define ip'^{x) — ip{x), then 

3 

^•"(x) = Vi*=-i(a;) -f E "™^™(^)^''"^(^ + ™) 

m— 1 

3 

+ amU^x - m)'ij/-\x - m) , (4) 

m—l 



and similarly for 'ip'^{x), with k = 1 . . . K. Here am is a 
strength factor which controls the amount of the smear- 
ing in spatial direction m. To keep the magnitude of the 
resulting correlation functions numerically under control 
it is useful to apply a rescaling factor after each smear- 
ing step. A suitable choice is motivated by the following 
observation: For the sake of this argument only, inter- 
pret ipj V' as Hilbcrt space operators obeying standard 
anti-commutation relations, then we find 

3 

{^\x)JHx)}^l{l + 2j2al), (5) 

m— 1 

where 1 is the color-spin unit matrix. Thus, with 

3 

i?=(l + 2E«m)^ (6) 
m—l 

the smeared quark fields are rescaled as 

after each iteration step Q. 

Similarly gauge field fuzzing is implemented defining 
Um{x) — Um{x), and then adding a sum over staples 

Ut{x) = Ut-\x) (8) 

3 

+ [Un-\x)Ut-\x + n)Ut^^(:x + m) 

n=l 

-rUt^^x - n)Ut\x - h)U^-\x + m-h)], 

with fuzzing strengths pm in directions m — 1,2,3. 
Again, for numerical reasons, it is useful to rescale the 
fuzzy link variables. These are no longer elements of 
SU(3). In particular, 

\Tr{u'^,\x)Ut{x)^ ^ S\ (9) 

which suggests the normalization prescription 

Ut{x) ^ S-'ut{x) (10) 

enforcing a trace of one after each fuzzing step 

Both normalization prescriptions are clearly invariant 
under gauge transformations. We use common strength 
parameters a™ = pm = 2.5 for all spatial directions. 

The basis of the current simulation are meson opera- 
tors Ox{k}{t) built from (P) and where k = 1. . .K 
indicates that the quark and gauge fields have been re- 
placed with smeared and fuzzed fields Q and JS)), respec- 
tively, at a common level k. The corresponding, single 
meson, correlation functions 

Cx{Ki}{iM= (11) 

{Oxiu}{t)0^x{i}ii^)) - {Ox{k}{t)){0^x{i}it^)) 

are the elements of a. K x K matrix. We use if = 3, 
and unsmeared fields {k — 0) are ignored. The same 
operators are used at both source and sink, so Cx(i,to) 
is hermitian. Due to our choice of fiavor structure l(2Jl, 
the separable terms in (|ll|l are zero. 
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III. ANALYSIS 

The time evolution of the eigenvalues of Cx{t,to) de- 
termines the mass spectrum for the meson X. Toward 
extracting the latter a standard procedure would be to 
solve a generalized eigenvalue problem 

Cx{t,to)^it) = Cx{ti,to)^it)A{t) (12) 

where ti is fixed, "^{t) is a K x K matrix, its columns 
being the generalized eigenvectors, and A{t) is real diago- 
nal. This approach is equivalent to a redefinition (linear 
combination) of the operators that make up the corre- 
lator matrix. Its purpose is to improve the numerical 
'quality' of effective mass function plateaus derived from 
the eigenvalues. 

We shall not rely on effective mass functions here, but 
rather analyze correlator eigenvalues by means of the 
maximum entropy method. Within this framework, solv- 
ing of the generalized eigenvalue problem H12(l offers no 
advantage from a numerical point of view. Thus we will 
diagonalize Cx{t,to) directly. 

Many eigenvalues, behaving exponentially with i, 
rapidly vanish into numerical noise. Conventional diag- 
onalization methods do not work well under those cir- 
cumstances. Singular value decomposition (SVD), on the 
other hand, is ideally suited to the task 16]. The SVD 
reads 

Cx{t,ta) ^ Ux{t,to)^x{t,to)vl{t,to), (13) 

where Ux{t, to) and Vx{t, to) are unitary in our case [25l| . 
and ^xit,to) = dia,g{ax;iit, to) ■■■ ax-Kit, tg)) contains 
the singular values satisfying <Tx;k{t,to) > 0. The rela- 
tion of the singular values (Jx;k{t,to) to the eigenvalues 
^X;k{t,to) of Cx{t, to) is given by the following theorem: 
If Cxit^to) is non-degenerate and positive semi-definite 
then the set of eigenvalues {Xx-k{t,to)\k = 1 . . . K} and 
singular values {crx-^kit,to)\k — I . . . K} are the same. 
Proving this is a simple exercise in linear algebra, the 
columns of Uxit, to) and Vx{t, to) are eigenvectors which 
may differ by phase factors only. In case of a degener- 
acy, which is highly improbable in a numerical setting, 
the usual ambiguities with respect to eigenvectors apply. 
For simplicity we will refer to ax-k{t, to) as eigenvalues. 

As an example, in Fig. ^ we show all eigenvalues of 
Cx{t,to) for the hybrid meson operator, X = h, com- 
puted at the smallest pion mass, i.e. the largest k value. 
Diagonalization is done independently on each time slice. 
Note that the eigenvalues are separated by almost three 
orders of magnitude. 

Bayesian analysis methods have been considered for 
lattice QCD [H [H ^ il HI HI , but are not widely 
used despite obvious advantages. For instance, one may 
maximize the time slice fitting range. The time slices 
used for all correlation functions in this work appear in 
Fig. n as filled plot symbols. 

The Bayesian fitting problem can be set up and solved 
in a variety of ways. We here adopt the framework de- 
scribed in detail in Ref. [23j |. A brief summary will be 
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FIG. 1: Eigenvalues of the 3x3 correlation matrix Cx(t,to) 
for the hybrid meson X = h a.t the smallest pion mass. Error 
bars are statistical, derived from a jackknife procedure [l^ . 
The lines are MEM fits as explained in the text. Open plot 
symbols indicate data points not included in the analysis. 

useful: The lattice provides correlation function data, say 
a{t), denoting any set (Jx;k{t,'to) on the given time slice 
range t = 2 ... 22. Our data should be well described by 
the model 

/"OO 

F(p|t) = / duo p{uj) cosh{uj{t ~ t^)) , (14) 
Jo 

where tc = 12 and p{uj) is a spectral density function. 
From the Bayesian perspective the numbers p{oj) are in- 
terpreted as stochastic variables. For these, employing 
Bayes' theorem, a conditional probability distribution 
function P[p <— a] is constructed. Aside from inessen- 
tial normalization it has the form 

P[p ^ a] « P[a ^ p\P[p\ . (15) 

Here P\a ^ p] is known as the likelihood function. We 
construct it from the x^-distance between the data and 
the model 

= E (^(^i) - ^(H^i)) r-\t,M) {^{h) - Fipm , 

(16) 

with T{ti,t2) being elements of the covariance matrix. 
Then P[a <— p] = exp(— x^/2) is the choice for the likeli- 
hood function. Bayesian inference uses information one 
might have about data, for example, a physical upper 
limit on the spectral masses. Information like this is 
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contained in P[p], known as the Bayesian prior. If no 
prior information is available, except say a physical mass 
range, a suitable choice is based on the Shannon-Jaynes 
entropy 



900 



S = I duj [ p{i^) — m{uj) — p{uj) In ^ 







m(uj) 



(17) 



Here J7 is the mentioned cutoff, and the function to(w) is 
called the default model. It provides a reference point for 
the spectral density in the sense that 5 < 0, with 5* — if 
p{lli) = m{ijj). Then P[p\ = exp(Q;S') is the choice for the 
Bayesian prior, introducing a new parameter a. Hence 
the conditional probability ifTCjl for the spectral density 
p becomes 



P[P 



(T oc e 



-(xV2-aS) 



(18) 



The idea then is to find a spectral density function p 
which maximizes P[p <— a], the posterior probability, at 
a fixed data set a. Bayesian inference within this frame- 
work is known as the maximum entropy method (MEM). 
In it was demonstrated that this problem can be 
solved by simulated annealing or cooling. Consider the 
partition function 

Zw= I [dp\e~P'^^^p'^ with W[p\ = (19) 



It involves the generation of equilibrium configurations 
[p] while gradually increasing /S^y from an initially small 
value, following some annealing schedule. For details we 
refer the reader to Ref. |23|. We here only mention that 
the resulting spectral density p(iJ) is extremely insensi- 
tive to both the default model m{uS) as well as to the 
entropy strength parameter a. 

The numerical implementation of the above scheme re- 
quires discretization of the lo integrals in ifT^ and ifTTjl . 
Choosing ajSl = 2.4 for the spectral mass cutoff in (|17|) 
we use otAcj — 0.05. The default model is a constant 
function m(a;) = lO^^at for all data sets crx;fe, and the 
entropy strength a is slightly adjusted, in each case, so as 
to render the ratio of aS to between 0.1 and 0.01, 
for the final p(io). Again, those parameter choices may 
be varied by several orders of magnitude without signif- 
icantly altering the resulting mass spectrum, see 23] for 
a discussion. 

The spectral density functions p(yi) for the eigenvalues 
fh;fc, k = 1,2,3, of the hybrid meson correlation matrix 
are displayed in Figs. EEEl respectively. The four panels 
in each of the figures show the results for the four hop- 
ping parameter values as listed in Tab. ^ below. The top 
panel, in each case, belongs to the largest k, or the small- 
est pion mass, and give rise to the model fits shown as 
solid lines in Fig.^ Each spectral density comes from 16 
random starts of the annealing process (|19|l . The center 
histogram (thick solid line) represents the average while 
the envelope histograms (thin dashed line) indicate the 
standard deviation with respect to the cooling starts. 
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FIG. 2: Spectral density functions (thick solid histogram 
lines) for the first (largest) eigenvalue set crji.fc, = 1, of the 
hybrid meson correlation matrix. The four panels correspond 
to increasing pion masses (decreasing hopping parameters k) 
from top to bottom. The envelop histograms (thin dashed 
lines) indicate errors as explained in the text. Smooth solid 
lines correspond to Gaussian fits. 



In most cases the spectra exhibit isolated peaks, say 
bn = {uj\u: S peak Then one may compute, for 

each peak, n = 1,2, the volume Z„, the mass E'n, and 
the width A„ , according to 



diopiyi) 



En. — Z„ 



du)p(iS)LL} 



(20) 
(21) 
(22) 
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FIG. 3: AsFig.|21 but for the second (intermediate) eigenvalue FIG. 4: As Fig. El but for the third (smaUest) eigenvalue set 
set ah,k, k = 2. ah,k, k ^ 3. 



Those (average) quantities constitute the information 
that can reasonably be expected to result from the 
Bayesian analysis. 

If overlapping peaks are present the above analysis 
gives inaccurate results. However, one may then resolve 
the spectral peaks by fitting a superposition of Gaussians 
to p{ol!)- Because we are forced to do this in at least one 
case, see Fig. |21 we have consistently made fits with one 
or two Gaussians, as required, of all spectral densities for 
the hybrid meson. To do this a Levenberg-Marquardt al- 
gorithm was employed, and the initial set of parameters 
were conveniently taken from applying H2()|I "H21 |I . The 
results of those fits are shown as smooth solid lines in 
Figs.EEIl 

A complete list of the results is given in Tab.HI There, 



we also show the MEM results for standard mesons 
X =^ TT, p, ai, which were obtained in the same manner, 
i.e. from 3x3 correlation matrices, etc. The masses of 
those mesons are used for reference to the physical scale, 
and extrapolation to the physical pion mass region, as 
discussed in the next section. The peak volumes Z„ are 
related to matrix elements {n\^k{t)\0), where ^k{t) is a 
linear combination of the (three) operators used to con- 
struct the correlator matrix . As such they are a mea- 
sure for how strongly the set of operators excite states |n) 
from the vacuum. The overall normalization, however, is 
arbitrary. 

Some comments on uncertainties are in order: As a 
matter of course jackknife configurations [T^ were gener- 
ated to estimate statistical errors. Those are reflected in 
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TABLE I: Results ol the MEM analysis of correlator matrix 
eigenvalues for meson operators X at four values of the hop- 
ping parameter k. The number k is the eigenvalue label. The 
entries for the peak volume Z„, the peak energy _E„, and the 
peak width An, for the primary peak n = 1 and secondary 
peak n — 2, where appropriate, are the results of Gaussian 
fits to the computed spectral density functions, see Figs. 
13 The uncertainties are standard deviations resulting from 
16 independent cooling starts. The presence of ±3 indicates 
factors of 10=^^ 
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K 


k 


n 




at En 




TV 


0.140 


1 


1 


18.4(0.8) 


0.50(1) 


0.0577(1) 




0.136 


1 


1 


17.0(1.4) 


0.60(1) 


0.0577(6) 




0.132 


1 


1 


16.1(1.4) 


0.71(1) 


0.057(1) 




0.128 


1 


1 


14.9(1.2) 


0.82(1) 


0.055(3) 


P 


0.140 


1 


1 


46.0(1.9) 


0.535(7) 


0.072(8) 




0.136 


1 


1 


47.1(0.7) 


0.628(8) 


0.067(6) 




0.132 


1 


1 


44.4(1.3) 


0.725(7) 


0.070(8) 




0.128 


1 


1 


42.6(0.5) 


0.828(4) 


0.067(5) 


ai 


0.140 


1 


1 


30.1(1.0) 


0.69(l)o 


0.13(1) 




0.136 


1 


1 


21.7(0.9) 


0.77(l)o 


0.11(1) 




0.132 


1 


1 


17.2(0.4) 


0.875(9)o 


0.105(8) 




0.128 


1 


1 


12.8(0.3) 


0.979(9)o 


0.096(9) 


h 


0.140 


1 


1 


0.10(6)+3 


0.62(4) 


0.16(1) 




0.136 


1 


1 


0.08(2)+3 


0.66(3) 


0.127(9) 




0.132 


1 


1 


0.07(l)+3 


0.76(2) 


0.079(4) 




0.128 


1 


1 


0.09(l)+3 


0.86(2) 


0.083(4) 


h 


0.140 


1 


2 


0.10(5)+3 


0.95(9) 


0.3(2) 




0.136 


1 


2 


0.123(9)+3 


1.02(2) 


0.17(2) 




0.132 


1 


2 


0.106(8)+3 


1.13(2) 


0.17(3) 




0.128 


1 


2 


0.079(6)+3 


1.26(2) 


0.21(3) 


h 


0.140 


2 


1 


0.455(8) 


0.756(7) 


0.26(2) 




0.136 


2 


1 


0.449(9) 


0.837(9) 


0.27(2) 




0.132 


2 


1 


0.421(5) 


0.900(6) 


0.26(1) 




0.128 


2 


1 


0.395(8) 


1.002(8) 


0.29(2) 


h 


0.140 


3 


1 


1.28(8)-3 


0.662(8) X 


0.11(3) 




0.136 


3 


1 


1.6(l)-3 


0.771(9)x 


0.09(3) 




0.132 


3 


1 


1.7(2)-3 


0.872(8) X 


0.09(3) 




0.128 


3 


1 


1.77(8)-3 


0.971(5) X 


0.09(2) 


h 


0.140 


3 


2 


6.3(2)-3 


1.395(8) 


0.26(4) 




0.136 


3 


2 


6.1(2)-3 


1.58(1) 


0.27(3) 




0.132 


3 


2 


5.9(l)-3 


1.72(1) 


0.32(2) 




0.128 


3 


2 


5.7(2)-3 


1.87(1) 


0.28(4) 



the error bars of the eigenvalue correlator plots of Fig.^ 
Hovifever, in the ensuing MEM analysis jackknife errors 
are hard to maintain, because extraction of the spectral 
density by way of simulated annealing is very time in- 
tensive. (Ideally, the cooling steps should be infinitely 
small.) For this reason the MEM analysis was only ap- 
plied to gauge configuration averages. Nevertheless, the 
peak widths of the spectral density functions naturally 
emerge as a measure of uncertainties for the masses. One 
may even argue that those are the 'true' measures of un- 
certainties if compared to statistical (jackknife) errors, 
because the latter can in principle be made arbitrarily 
small by increasing the number of gauge configurations, 
even though the lattice data may be poor. Peak width 
uncertainties and plateau statistical errors are hard to re- 



late, one reason being that they stem from different data 
sets. This notwithstanding we observe that the masses 
of Ref. 0, which carry statistical errors, are generally 
consistent with the ground state meson masses obtained 
in the current analysis using peak width uncertainties. 

IV. DISCUSSION OF RESULTS 

In order to set the spectral masses into a physical con- 
text we have attempted extrapolations to the physical 
pion mass region. Thus we consider the computed p, ai 
and h spectral masses as functions of . To our knowl- 
edge, predictions for this dependence from chiral pertur- 
bation are not available, as a hybrid meson is in- 
volved. We therefore employ a heuristic model. It has 
been used in Ref. 0, where also some motivation was 
given. For a spectral mass M = atm we attempt fits 
with the model 

M = p + qx + r ln(l + x) with x = {atm-,^)'^ , (23) 

using three parameters p, q, r. The linear term is stan- 
dard in xPT inspired models while the logarithmic term 
is purely heuristical but happens to yield excellent fits to 
our data. The results are illustrated in Fig. |S1 

The extrapolated p meson mass, at x — 0, is used to de- 
termine the physical scale for this simulation. We obtain 
at = 0.36(6) GeV"^ = 0.072(11) fm {rup = 775.8 Me V). 
Tableimthen gives a list of the extrapolated masses at Ex 
along with the physical values Ex- Also listed are esti- 
mates of the peak widths at x = 0. Those are obtained 
by randomizing the data, replacing the four data points 
with at En + ^at A„ where ^ is a normal distributed ran- 
dom deviate with variance one, and then repeating the 
fits with the model H23|l obtaining random values for the 
extrapolated spectral masses atEx- A few thousand ran- 
domizations then give the results Ax in Tab. ^1 Uncer- 
tainties on Ax, being second order, were not calculated. 
The errors given in parentheses are obtained in the same 
way, but using the uncertainties from Tab. instead of 
the A„. Unfortunately, extrapolation tends to amplify 
the peak widths because they become larger at smaller 
pion masses. As a measure of uncertainty the extrapo- 
lated Ax, in Tab.im thus seem less useful. 

As a first observation we note that the ai meson mass 
in Tab.^]comes out very close to the experimental value 
1230 MeV, which is the mass of the ai(1260) meson, in 
the nomenclature of 0. We take this as a hint that 
the extrapolation model is adequate, and therefore, the 
extrapolated results for the h mass spectrum, in Tab.HTl 
may not be without merit. 

In the light of this we continue to observe that the 
energy levels of Xq^^^) for and ^(2,1) in Tab. ^ 

match the experimental masses 1376 MeV and 1653 MeV 
of the 1 resonances known as 7ri(1400) and 7ri(1600) 
in Ref. 2]. One should not overinterpret the closeness of 
the computed masses to the experimental levels, because 
systematic errors associated with the extrapolation could 
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FIG. 5: Plots of spectral meson masses versus the squared 
pion mass x = (atnin)^ and fits with the model (1231 . Solid 
lines refer to masses from primary (n = 1) spectral peaks, 
and dashed lines to secondary (n = 2) peaks, see Tab.Q] The 
extrapolation of the p meson mass to a; = is used to set the 
physical scale. 



note that the extrapolated spectral level sequence only 
emerges beyond the level crossing of the states and 
ai at around x ~ 0.1, see Fig. 

Another remarkable results concerns the level /i(3,i). 
Its energy spectrum coincides with the one generated by 
the ai meson operator. In Tab. ^ the symbols o and 
X are attached to the corresponding levels. These are 
also used in Fig. [S] where an almost perfect match is 
apparent, only the data point at the smallest pion mass is 
slightly off. It causes, however, unrealistic extrapolation 
results, see the entries in Tab. |n] marked with a x. The 
key to an explanation of this curious result goes back to 
the observation made in Sect, m namely that parity is 
guaranteed only in the limit of a large number of gauge 
configurations. The level /i(3,i) comes from a correlator 
eigenvalue about four orders of magnitude less than the 
dominant one. It may thus single out the contamination 
of the lattice signal with a 1+"'" state, which happen to be 
the quantum numbers of the oi meson Q]. Although this 
observation is not directly relevant to our current interest 
in 1 exotic states, it nevertheless adds credence to the 
current simulation and analysis method. 

Lastly, the largest levels /i(i,2) and /i(3,2) in Tab. ^ 
point at a mass somewhat above 2 GeV. We speculate 
that at least one of those levels coincides with the 1 ^ 
resonance at 1.9 GeV uncovered in Ref. 0. There the 
space of operators was larger, including a Trai two-hadron 
field in addition to h, and this could be the cause for a 
lowering of the energy level from 2 GeV to 1.9 GeV. If this 
interpretation is correct, then this state likely is a mix- 
ture of a hybrid meson and a two-meson state. We take 
this result as supplementary evidence to the outcome of 
Ref. 0. 



TABLE II: Extrapolated spectral masses Ex and peak 
widths Ax, for mesons X. The eigenvalue label is k and 
the spectral peak number is n, as in Tab. 13 The values for 
Ax and for the uncertainties (in parentheses) are obtained by 
randomization of the data points, as explained in the text. 



X 


k 


n 


at Ex 


atAx 


Ex[GeY\ 


Ax [GeV] 


p 


1 


1 


0.28(04) 


0.08 


0.7785 


0.22 




1 


1 


0.45(06) 


0.17 


1.23(0.17) 


0.47 


h 


1 


1 


0.52(19) 


0.37 


1.43(0.53) 


1.03 


h 


1 


2 


0.79(37) 


0.65 


2.19(1.03) 


1.80 


h 


2 


1 


0.59(04) 


0.27 


1.63(0.12) 


0.76 


h 


3 


1 


0.34(05) X 


0.20 


0.94(0.13) X 


0.55 


h 


3 


2 


0.84(07) 


0.11 


2.32(0.18) 


0.30 



be large. On the other hand the spectral level pattern of 
two 1 ^ states below 2 GeV is a unequivocal result of 
this simulation, which in no small part relies on the MEM 
analysis of the correlator data. It is also interesting to 



V. CONCLUSION 

Using the maximum entropy method, five distinct 
spectral levels have been uncovered for the J^*^ = 1 
exotic meson. Two of the spectral levels correspond with 
the 7ri(1400) and 7ri(1600) from 0. Two more levels 
possibly correspond with a resonance energy of 1.9 GeV 
previously determined by a decay width calculation us- 
ing Liischer's method 0. A fifth spectral level, at higher 
pion masses, tracks consistent with an operator repre- 
senting the ai(1260) meson, and we take this level to be 
a consequence of inexact parity symmetry. 

All of these spectral levels rely on extrapolations to 
TOtt — from relatively heavy pions. Although this may 
give rise to large systematic errors, the fact that the oi 
extrapolation came very close to its experimental value 
leads us to conclude at least two spectral levels for the 
1 ^ exotic meson will be below 2 GeV. 
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